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, Abstract: A study has been made of the physics interest and feasibility of ex- 

periments with nuclear beams in HERA. It is shown that such experiments widen 
considerably the horizon for probing QCD compared to that from free nucleon tar- 
gets. In addition there is some sensitivity to physics beyond the standard model. 
Hence the option to include circulating nuclear beams in HERA allows a wide range 
\Q . of physics processes to be studied and understood. 
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^ ■ 1 Introduction 

> : 

The successes of QCD in describing inclusive perturbative phenomena have moved the focus of 
■ investigations to new frontiers. Three fundamental questions to be resolved are the space-time 
structure of high-energy strong interactions, the QCD dynamics in the nonlinear, small coupling 
domain and the QCD dynamics of interactions of fast, compact colour singlet systems. 

The study of electron-nucleus scattering at HERA allows a new regime to be probed exper- 
imentally for the first time. This is the regime in which the virtual photon interacts coherently 
with all the nucleons at a given impact parameter. In the rest frame of the nucleus this can be 
visualized in terms of the propagation of a small qq pair in high density gluon fields over much 
larger distances than is possible with free nucleons. In the Breit frame it corresponds to the 
fact that small x partons cannot be localized longitudinally to better than the size of the nu- 
cleus. Thus low x partons from different nucleons overlap spatially creating much larger parton 
densities than in the free nucleon case. This leads to a large amplification of the nonlinear 
effects expected in QCD at small x. The HERA ep data have confirmed the rapid increase 
of the parton densities in the small x limit predicted by perturbative QCD. However the lim- 
ited x range available at HERA makes it difficult to distinguish between the predictions of the 
DGLAP evolution equations and the BFKL-type dynamics. Moreover, the nonlinear effects 
expected at small x are relatively small in ep scattering in the HERA kinematic domain and 
it may be necessary to reduce x by at least one order of magnitude to observe unambiguously 
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such effects. However, the amplification obtained with heavy nuclear targets allows an effec- 
tive reduction of about two orders of magnitude in x making it feasible to explore such 
nonlinear effects at the energies available at HERA. The question of nonlinear effects is one 
of the most fundamental in QCD. It is crucial for understanding the kind of dynamics which 
would slow down and eventually stop the rapid growth of the cross section (or the structure 
function, F 2 ) at small x. It is also essential in order to understand down to what values of x the 
decomposition of the cross section into terms with different powers of remains effective. It 
is important for the understanding of the relationship between hard and soft physics. One can 
also study the dynamics of QCD at high densities and at zero temperatures raising questions 
complementary to those addressed in the search for a quark-gluon plasma in high-energy heavy 
ion collisions. 

Deep inelastic scattering from nuclei provides also a number of ways to probe the dynamics 
of high-energy interactions of small colour singlet systems. This issue started from the 
work of Gribov |IJ who demonstrated the following paradox. If one makes the natural (in 
soft physics) assumption that at high energies any hadron interacts with a heavy nucleus with 
cross section 2ttR 2 a (corresponding to interaction with a black body), Bjorken scaling at small 
x is grossly violated - o^a oc InQ 2 instead of To preserve scaling, Bjorken suggested, 
using parton model arguments, that only configurations with small p t < Pto are involved in the 
interaction (the Aligned Jet Model) M. However, in perturbative QCD Bjorken's assumption 
does not hold - large pt configurations interact with finite though small cross sections (colour 
screening) , which however increase rapidly with incident energy due to the increase of the gluon 
density with decreasing x. Hence again one is faced with a fundamental question which can 
only be answered experimentally: Can small colour singlets interact with hadrons with cross 
sections comparable to that of normal hadrons? At HERA one can both establish the x, Q 2 
range where the cross section of small colour singlets is small - colour transparency, and look 
for the onset of the new regime of large cross sections, perturbative colour opacity. 

Another fundamental question to be addressed is the propagation of quarks through 
nuclear matter. At large energies perturbative QCD leads to the analogue of the Landau- 
Migdal-Pomeranchuk effect in quantum electrodynamics. In particular Baier et al. || find 
a highly nontrivial dependence of the energy loss on the distance, L, travelled by a parton 
in a nuclear medium: the loss instead of being oc L is oc L 2 . Several manifestations of this 
phenomenon can be studied at HERA. 

There is also an important connection to heavy ion physics. Study of eA scattering at 
HERA would be important for the analysis of heavy ion collisions at the LHC and RHIC. 
Measurements of gluon shadowing at small x are necessary for a reliable interpretation of the 
high p t jet rates at the LHC. In addition, the study of parton propagation in nuclear media is 
important for the analysis of jet quenching phenomena, which may be one of the most direct 
global signals of the formation of a quark-gluon plasma. 

Current fixed target data on lepton-nucleus scattering only touch the surface of all these 
effects due to the limited Q 2 range of the data at small x. Indeed the Q 2 range of these data 
is too small to distinguish the contribution of the vector meson dominance behaviour of the 
photon from its hard QCD behaviour at small x. The range of x and Q 2 in experiments with 
nuclei at HERA compared to the fixed target experiments is shown in Fig. [I]. It can be seen that 
at HERA the kinematic range will be extended well into the deep inelastic scattering region. 

To address the questions discussed above we identify the primary experimental programme 
for nuclei in HERA as: 



2 




10 " 5 10 " 4 10 " 3 10 " 2 10 



Figure 1: The kinematic region covered by experiments at HERA compared to fixed target 
data (shaded region). 

• Study of the x and Q 2 dependence of nuclear shadowing over a wide Q 2 range. 
This will allow the processes limiting the growth of Fi as x tends to zero to be studied in detail. 

• To establish the difference between the gluon distributions of bound and free 
nucleons. This will allow the part played by gluon fusion in the shadowing process to be 
studied directly. 

• Study of diffractive processes: to see if the pomeron generated by nuclei shows any 
difference from that generated by free nucleons. Processes such as vector meson production 
can also be used to search for colour transparency. 

• Study of hadronic final states. This allows the propagation of partons in the nuclear 
medium to be studied as well as the multiplicity fluctuations discussed later. 

The proceedings of Working Group 8 are organised as follows. First we give an experimental 
overview in which we demonstrate the feasibility of carrying out this experimental programme. 
Then we give a theoretical overview in which we explain the relevance of the programme to 
QCD. Finally, we give the detailed contributions on different topics which demonstrate the 
depth of the physics interest. The proposed measurements in the main will be possible with 
the existing detectors HI and ZEUS measuring down to low Q 2 and with luminosities at the 
level of 1-10 pb _1 per nucleon. The contribution of Chwastowski and Krasny |4]] shows that if 
the detectors could extend their rapidity coverages various experiments of interest to nuclear 
structure physicists become feasible. 



2 Experimental Overview 
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2.1 Introduction 



In the following subsections the feasibility of the measurements defined above as the primary 
experimental programme is investigated. The nuclear targets should each have Z/A of 1/2. 
Hence the energy of each nucleon in a deep inelastic collision will be half that of the HERA 
proton energy of 820 GeV i.e. 410 GeV. The electron energy is assumed to have the standard 
value of 27.6 GeV. We show that most studies can be carried out with the existing detectors 
requiring luminosities between 1-10 pb" 1 per nucleon. While the possibility of storing heavy 
nuclei up to Sn and Pb is very attractive, a program covering the light isoscalar nuclei (D, 4 He, 
C, S) would by itself have a major discovery potential. The necessary radiative corrections are 
described in the contributions of Kurek || and of Akushevich and Spiesberger || who show 
that such corrections can be kept under control using suitable cuts on the data. In most of the 
experiments it is proposed to measure ratios of yields. Hence, in order to minimise systematic 
errors, it is desirable to store different nuclei in HERA simultaneously. Otherwise frequent 
changes of the stored nucleus in the beam will be necessary. 



2.2 Shadowing Measurements Using Nuclei in HERA 

The accuracy of shadowing measurements for an experiment in which nuclear targets are stored 
in HERA is calculated. It is shown that for luminosities of 2 pb -1 per nucleon the measurements 
would extend considerably the accuracy and range of the existing data. 

The x dependence of the differences in the nucleon structure function between bound and 
free nucleons has been well measured in recent years in fixed target experiments [0, |8|, |9], [K], 
pTTfl following the discovery of the differences in the 1980s [12, |13|, [14], [15], 16| . The present 
experimental situation on this x dependence is briefly summarised in Fig. 0. However, the Q 2 
dependence is not well measured. In the shadowing region the data are at such low Q 2 values 
that they are arguably not even in the deep inelastic regime. 

There are many different models |T7|] for the effects in the different regions shown in Fig. ^] 
which are all compatible with the existing data. In the shadowing region the Q 2 range of the 
data is insufficient to separate the different contributions from the vector dominance behaviour 
of the photon and QCD effects such as parton fusion. Measurements over the extended x and 
Q 2 ranges, which would become possible at HERA, will give more information to help separate 
the models and help us understand the phenomena which limit the rise of the nucleon structure 
function F 2 at small x; e.g. see references []18" , |i~9 . 



Is it feasible to study shadowing in HERA? 

To answer this question we assume that nuclear beams can be stored in HERA and lumi- 
nosities of 2 pb -1 per nucleon, shared between 2 nuclear targets, can be achieved (i.e. 2/ A pb -1 
per nucleus, where A is the atomic weight). With this definition of luminosity, rate computa- 
tions should use cross sections per nucleon. The counting rates in bins of Q 2 and x are then 
estimated to assess the statistical errors on the measurements of the ratios of F 2 A / F 2 D where 
the nuclear targets are assumed to be He, C or S. The cross sections are calculated from the 
one photon exchange formula 
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Figure 2: Partial compilation of results for F 2 / F® (from |]]). Note that below x ~ 0.01, the 
average Q 2 value of the data is smaller than 2 GeV 2 . 



The nucleon structure functions, F 2 and R, are computed from the MRS (A) set |2(J. Deep 
inelastic scattering events are assumed to be detectable with 100% efficiency if the scattered 
electron energy E' > 5 GeV and its scattering angle is more than 3 degrees to the electron 
beam. Alternatively events are assumed to be detectable with 100% efficiency if a quark is 
scattered out by an angle of more than 10 degrees to the proton beam and with an energy more 
than 5 GeV. Nuclear effects on the structure functions are neglected. Such effects, which are 
at the 10% level, will have little effect on the statistical accuracy of measurements of ratios of 
structure functions, F 2 A /F 2 D . 

Fig. [3] shows the estimated statistical errors on the ratios (averaged over Q 2 ), with these 
assumptions, together with the measurements of the NMC |§. In many cases the statistical 
errors are < 1% i.e. smaller than the sizes of the points. It can be seen from Fig. [3] that this 
statistical precision will allow high accuracy measurements of the shadowing ratios down to 
lower x values than in fixed target experiments and over a much wider range of Q 2 . 

Fig. [| shows the statistical precision of the slopes d(F^/F 2 D )/dlnQ 2 estimated at HERA 
compared to the NMC data. Impressive precision is possible at HERA, presumably due to the 
much larger Q 2 range covered. Measurements over such a large Q 2 range will allow the precise 
predictions of the parton fusion model to be tested. 

If at least two different nuclear targets can be stored in different bunches in HERA during 
experimental running the systematic errors should be similar in magnitude to those in fixed 
target experiments JFj. Radiative corrections for deep inelastic scattering from heavy nuclear 
targets will be necessary and these will be applied with the appropriate cuts on the data as 
discussed in these preceedings || ||. 
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In conclusion, measurements of the ratios (F^/F^) at HERA will extend the data to smaller 
values of x and much larger values of Q 2 than in fixed target experiments. Impressive precision 
on the Q 2 dependence will be possible from which the mechanism which leads to the limitation 
in the rise of F2 at small x and large Q 2 can be studied. 
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Figure 3: Ratio of the nucleon structure function in carbon to that in deuterium as a function 
of x. The NMC data |§ (open squares) are shown in comparison to data with the estimated 
statistical accuracy of an experiment of luminosity 1 pb~ x per nucleon at HERA. 



2.3 The Accuracy of The Gluon Density Measurements Using Nu- 
clei in HERA 

2.3.1 Introduction 

The major contribution to shadowing from nonlinear QCD effects is thought to arise from 
multigluon interactions such as gluon fusion effects. Such effects are amplified at higher x in 
nuclei due to the larger target size so that they should become visible in the HERA kinematic 
range. Similar effects are expected to limit the growth of the nucleon structure function F 2 
at high Q 2 and low x. Hence it is interesting to look for such effects directly on the gluon 
distribution by looking for differences between the gluon density in bound and free nucleons. 
Different ways of doing this are studied below. 
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Figure 4: The slopes d(F% / F®) / d\nQ 2 as function of x showing the NMC data [[jj] (open 
squares) and the statistical accuracy of an experiment with 1 pb _1 per nucleon at HERA. 






Figure 5: The processes giving rise to 2 + 1 jet topology. P denotes here the projectile particle, 
a proton or a nucleus, Q 2 is the four-momentum transfer and x is the Bjorken variable. 



2.3.2 Determination of the Gluon Density in the Nucleus from the Jet Rates 

In the majority of large Q 2 deep inelastic electron-nucleus scattering events, at HERA, the 
scattered quark and the remnant of the nucleus will hadronize to form two jets. Such a topology 
is often called the 1 + 1 jet configuration. Occasionally, however, more jets will be produced. 
In Fig. |5] the partonic processes giving rise to the 2 + 1 jet topology are shown. These processes 
are called "boson-gluon fusion" (process a) and "the QCD Compton scattering" (processes b). 

The 2 + 1 jet events can be related to the partonic processes shown in Fig. || if the invariant 
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mass of the system of two jets s = (pjeti + Vjet2) 2 is significantly larger than the typical scale 
of the strong interactions, so that perturbative QCD can be used. The contribution of these 
processes to the total cross section depends upon the value of the coupling constant ot s defining 
the strength of quark-gluon coupling and upon the momentum distribution of the incoming 
gluon (quark). The fractions of the parent bound nucleon momentum carried by the incoming 
partons, x Pt9 , are constrained by the value of x, the total hadronic mass, W, and the invariant 
mass of the two jet system s: 

x Pt9 = x + s/W 2 . (2) 
For values of s/W 2 > 0.01, the values of x p>g must be large. In this kinematic domain the par- 



tonic distributions in bound nucleons have been well measured in fixed target experiments []T7 
Thus, the coupling constant, a s , can be derived from the measured rate of 2 + 1 jet events. 
In turn, we shall be able to use this a s value to determine the gluon momentum distribution 
at smaller x g , corresponding to small values of both x and s/W 2 (note that at small x the 
contribution of the QCD Compton processes (Fig. ||b,c) to the total "2+1" jet cross section is 
expected to be small). Such an analysis has been recently carried out by the HI collaboration 
f2~Tf using deep inelastic electron-proton scattering data collected at HERA. The systematic er- 



rors of the resulting gluon density are large and to some extent uncertain. They are dominated 
by the uncertainties in relating the measured rate of observed jet topologies to the basic QCD 
processes involving quarks and gluons. Unfolding the gluon densities involves modelling the 
hadronisation of the quarks which is necessary for simulating the detector effects but is only 
weakly constrained by the data. In addition several jet algorithms can be used leading to differ- 
ences in jet counting. There exist as well ambiguities in the QCD calculation of the processes 
shown in Fig. ||. The amplitudes of the processes 7*g — > jeti + jet 2 and 7*g — > jeti + jet2 
have poles corresponding to the collinear emission of jets with respect to the direction of the 
incoming gluon (7*) and to the emission of jets of small invariant mass. Since one must use 
fixed order perturbative QCD, jets reconstructed in phase space close to the poles must be 
avoided to diminish the sensitivity to higher order terms which have so far not been calculated. 
This can be achieved by using only jets of high invariant mass (s > 100 GeV 2 ) and by rejecting 
events in which one of the jets is emitted at a small angle with respect to the incoming proton 
direction. In the region of small s the jet rates calculated using leading order and next to lead- 
ing order approximations are significantly different [^] indicating that higher order corrections 
are necessary for an unambigous determination of the gluon density. 

Will it be possible to use the jet method to determine the gluon density in nuclei given 
the uncertainties described above and at the same time adding the extra ambiguity related 
to jet formation in the nuclear medium? Most likely it will be difficult to obtain satisfactory 
precision in measuring the absolute gluon distributions in the nuclei. However, we expect that 
good precision can be achieved in measuring the ratios of gluon distributions for various nuclei. 
The uncertainties due to the jet finding algorithms and to the modelling of the hadronisation 
processes will largely cancel in the ratios if high energy jets are used. High energy jets are 
expected to be formed outside the nucleus. In addition we shall be able to select jets in the 
restricted phase space region where the effects of rescattering of slow particles belonging to the 
jet but formed inside the nucleus are small. The energy loss of quarks and gluons traversing 
the nuclear medium prior to hadronisation is expected to be below 350 MeV/fm according to 
the estimation of [[| and should not give rise to large errors on the gluon density ratio. The 
expected statistical precision of the measurements of the ratio of the gluon densities for two 
nuclei of atomic numbers Al and A2 is shown in Fig. || It corresponds to a luminosity of 10 
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Figure 6: Statistical precision of the ratio of gluon densities in two 

xnG(xn, Ai)/xG(xn, A 2 ) at Q 2 = 30 GeV 2 determined from the "2 + 1" jet sample. 



nuclei 



pb~ l /A for each nucleus. In estimating this statistical precision we have followed the jet and 
kinematic region selection of and neglected all nuclear effects. We also assume that the 
contribution of the QCD Compton process can be unambigously subtracted. The x^ variable 
is the fraction of the bound nucleon momentum carried by the gluon. It is clear that good 
statistical precision on the measurement of the ratio of the gluon densities can be achieved at 
modest luminosities. 



2.3.3 The Gluon Distribution in Nuclei from Scaling Violations 

The gluon distribution can also be determined from the deviations from scaling of the structure 
function F 2 . We estimate the accuracy of such a determination of the gluon density for bound 
nucleons in an experiment with nuclear beams in HERA. The scaling violations of F 2 are 
strongly related to the gluon distribution of the target nucleon at small x values. The quantity 



^ ~ Iff w 

dlnQ 2 

is sensitive to the differences of the gluon distribution in bound and free nucleons and is roughly 
proportional to 5G/G where G is the gluon density at a particular x value and SG is the 
difference in gluon densities between bound and free nucleons. Hence it would be interesting 
to measure this quantity in an experiment with nuclear beams stored in HERA. 

To estimate the feasibility of such a measurement the accuracies of the determinations of 
the slopes dF 2 /d In Q 2 have been estimated from the expected counting rates for a 2 pb -1 
per nucleon run using the MRS (A) set of structure functions [P0[] . The slopes were obtained 
from a linear least squares fit to values of F 2 using the statistical errors calculated for such an 
experiment with the assumptions described in section |2.2j . The error in the quantity ip was 
then obtained assuming equal statistical errors for the nuclear and deuteron targets with a 1 
pb -1 per nucleon run for each. The error on the slope dF 2 / d\nQ 2 for the proton was neglected 
since this should be well determined from high luminosity proton running. 
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Figure 7: Accuracy of the difference between the slopes of F 2 between nuclei and deuterium 
as a function of x. The smooth curve shows the calculated slope for the free proton. 

The results are shown in Fig. [?]. For the purpose of this figure we assume that the gluon 
densities for bound and free nucleons are the same so that ^=0. The estimated statistical errors 
are then superimposed on these values and the smooth curve shows dF^/d In Q 2 for comparison 
computed from the MRS (A) set. Comparison of the errors on ip with the values of the slope 
for the proton shows that this method will be sensitive to differences in the gluon densities 
between bound and free nucleons of the order of 5 per cent of the total gluon density at x in 
the vicinity of 10~ 3 . Hence the measurement will be quite sensitive to nuclear effects on the 
gluon density. 
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2.3.4 Determination of the Gluon Density in the Nucleus from Inelastic J/ip pro- 
duction 

The expected accuracy of the determination of the bound to free nucleon gluon density ratio 
xG\a/xG\d using nuclei in HERA is estimated for inelastic photoproduction of J/ip mesons at 
Q 2 < 4 GeV 2 (eA -> eXJ/^j). 

Inelastic production of J/ip mesons has been used for a long time to extract or constrain 
the gluon distribution in the nucleon (see e.g. p3"|-P8|), assuming that the dominant mecha- 
nism is the photon-gluon fusion f2"9fl- [|3^| . Within this framework, the cross section is directly 



proportional to the gluon density. These calculations have been affected in the past by large 
normalisation uncertainties (up to factors of 2-5), which however have been greatly reduced 
recently fl33" |. Inelastic J/ip production dominates at values of z, the fraction of the photon 



energy carried by the meson in the nucleon rest frame, smaller than « 0.9. In the framework of 
the colour singlet model |32| , the gluon distribution is probed at a value of x, the fraction of the 
proton's momentum carried by the gluon, x = [m 2 ^/ z + pf /z(l — z)]/W 2 , where p t is the J/ip 
transverse momentum with respect to the virtual photon direction and W is the photon-nucleon 
centre-of-mass energy. The scale probed by this process is approximately m 2 ,^ ~ 10 GeV 2 . 
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Fig. | shows the expected statistical accuracy as a function of log 10 x for an integrated 
luminosity of 10 pb^/A Decays into e + e~ or [i + pT pairs have been assumed. The plot refers 
to the kinematic region Q 2 < 4 GeV 2 , z < 0.9. Nuclear effects have been neglected in the 
evaluation of statistical accuracies. 

Possible sources of systematic uncertainties are the luminosity, the branching ratio, the 
global acceptance (including trigger and reconstruction efficiency, muon or electron identifica- 
tion etc.), the feed-in from ip' production and the contamination from resolved photon events. 



In the recent HI |26| and ZEUS |2y| analyses, the total systematic uncertainty is approximately 
10-20%. By and large all the above contributions would cancel in a ratio for simultaneously 
stored nuclei, with the partial exception of the luminosity. In practice, for an integrated lumi- 
nosity of 10 pb -1 , the statistical uncertainty dominates. 
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Figure 8: Expected statistical accuracy (vertical error bars) on the ratio xG\a/xG\d as a 
function of log 10 x using inelastic photoproduction (Q 2 < 4 GeV 2 ) of J/ip mesons. An integrated 
luminosity of 10 pb~ 1 /A for each nucleus has been assumed. The horizontal bars indicate the 
size of the bins. 



2.4 Diffraction from Nuclei in HERA 
2.4.1 Introduction 



The measurements of the ZEUS and HI collaborations |34|, |35|] of deep-inelastic electron-proton 
scattering have revealed the existence of a distinct class of events in which there is no hadronic 
energy flow in an interval of pseudo-rapidity, rj, adjacent to the proton beam direction i.e. 
events with a large rapidity gap. Such events are interpreted as deep inelastic scattering from 
the pomeron, IP. Studies of events with a large rapidity gap from nuclear targets will allow the 
structure of the pomeron from a different source than the free nucleon to be determined. It will 
be interesting to see if these structures differ. In addition, the study of diffractive vector meson 
production will be interesting to search for the phenomenon of colour transparency. Such a 
phenomenon has not yet been convincingly seen although it is predicted in QCD. 



2.4.2 Expected Accuracy of the Measurements of the Pomeron Structure 

We commence by describing the terminology surrounding measurements of the Pomeron struc- 
ture. In our studies we shall use the four variables Pa, Q 2 , xa and tA, or equivalently Pa, Q 2 , 
xjp^a and tA, which are defined as follows: 

XA = 2p— q ' XF > A= q .p ; 0=-*; P A = 2q .(p-p>y ^ = ( p - p )-( 4 ) 
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Here q, P and P' are, as indicated in Fig. [|, the 4-momenta of the virtual boson, incident 
nucleus and the final state colourless remnant Y respectively. The latter can be either a 
coherently recoiling nucleus or any incoherent excitation of the nucleus carrying its quantum 
numbers. The variables in equation (f|) are related to each other via the expression: 

x A = p A X]p A . (5) 

We also introduce, in order to allow the comparison of measurements with nuclei of different 
atomic number A, the variables: 

x = x A ■ A xjp = x P!A -A (3 = (3 A t = t A . (6) 

Note that relation ([5]) rewritten in terms of the above variables still holds i.e. 

x = (3x F . (7) 

The A-rescaled variables can be directly related to the variables defined in [BBl B7J retaining 



their interpretation, as was given there, for processes in which only one nucleon of the nucleus 
interacts with the electron and the nucleon's Fermi momentum can be neglected. 

The variables x^a, %ip and (3 can be expressed in terms of the invariant mass of the hadronic 
system X, M x , the nucleus mass, Ma, and the total hadronic invariant mass W as 

Q 2 + Ml-t Q 2 + M 2 X fk . 
XF > A = Q 2 + W 2 -M 2 A W ' (8) 

Q 2 + M 2 X 

xp « ■ x, (9) 

P = Pa= Q*+ Q M x -t ^Q^W (10) 

In the kinematic domain which we shall consider here [M\ ^ ^ 2 an d Q 2 ■> \^\ < ^ M x ) 
xp,A m ay be interpreted as the fraction of the 4-momentum of the nucleus carried by the IP 
and (3 as the fraction of the 4-momentum of the IP carried by the quark interacting with the 
virtual boson. Note that in the interactions in which only one nucleon takes part, xjp can be 
interpreted as the fraction of the 4-momentum of this nucleon carried by the IP. 



We shall follow in our analysis the formalism defined in [gg, ^TJ an d introduce the diffractive 
structure function F^a as a function of three kinematic variables, derived from the structure 

function F^a- This latter structure function depends upon four kinematic variables chosen 
here as: x, Q 2 , xjp and t and is defined by analogy with the decomposition of the unpolarised 
total eA cross section. The total cross section can be expressed in terms of two structure 



functions F^f 1 and — : 2 ' A dia\ in the form 

" A ' \ 1 -y + ZZ „nu),„ tt ( F 2 X [f3,Q , Xjp, t), (11) 



D(4) 



dxdQ 2 dxjpdt xQ* { 2[l + J Rj (4) (/3,Q 2 ,x JP ,t)] 

in which y = Q 2 /sxa and s is the eA collision centre of mass (CM) energy squared. We shall 
discuss here the low y region and neglect the term containing pJa^' '■ 
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Figure 9: The diagram of the process with a rapidity gap between the system X and Y. The 
projectile nucleus is denoted here as p. 



We shall consider in the following measurements of Z A ^f — -, from which the structure 



function F^\ft, Q 2 , xp) = J F®^\f3, Q 2 , xjp, t) dt can be derived. The integration is over the 
range | t min \ <\t\<\t \u m where t min is a function of Q 2 , W 2 , ft and the mass of the the system 
Y, and \ t\u m is specified by the requirement that all particles belonging to the system Y remain 
undetected. The structure function will thus be derived from 

dV^eXY _ , 4vra 2 j y 2 \ n (3) 



The structure function F^\ft,Q 2 ,xp) can be measured at HERA using a sample of "ra- 
pidity gap" events i.e. events in which there is no hadronic energy flow over a large 7] interval. 
These events originate from coherent diffractive scattering {eA — > eA + Xdiff) and from inco- 
herent diffractive scattering (eA — > e(A — N) +N + Xdiff)- N denotes here the nucleons ejected 
from the incoming nucleus. The measurement of F^a {ft, Q 2 , xp) for several nuclei, separately 
for coherent and incoherent processes will provide an important test of the quark parton 
interpretation of diffractive processes and a unique means to find out how universal is the 
concept of the pomeron. 

We propose to measure the ratio of the structure functions 

RaiMP'&'Xf) = f ?,aK^ Q 2 , x p )/f£HH(P, Q 2 , Xjp), (13) 

where Al and A2 denote the atomic numbers of the two nuclei. This ratio can be measured at 
HERA with a very high systematic accuracy. The statistical precision of such a measurement 
will be of the order of 5 %, if luminosities of 10/A1 and 10/A2 pb™ 1 are collected for each 
nucleus. This is illustrated in Fig. ITUI where we show, as an example, Rai,az a ^ Q 2 = 12 GeV 2 



as a function of ft and x P . In order to estimate the statistical precision of the measurement we 
have used the RAPGAP |JS] Monte-Carlo and have assumed the event selection procedure of 
36], |3?|. We have not tried to model the nuclear dependence of the F^^ix, Q 2 , xjp) resulting 



in the ratio shown in Fig. ^ to be equal 1. 

Several distinct hypotheses concerning the deep inelastic structure of the diffractive pro- 
cesses can be verified (rejected) by measuring the Rai,A2'- 

• universal (independent of the source) pomeron structure and an A-independent pomeron 
flux leading to 

RaiMP'Q 2 i x ip) = i; (14) 
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Figure 10: The ratio Rai,A2(/3,Q 2 ,xip) plotted as a function of xjp for fixed values of (3 and 
Q 2 . The error bars correspond to a luminosity of 10/^41 and 10/A2 pb _1 for each nucleus. 



universal pomeron structure and an A-dependent pomeron flux leading to 

RaiM/3,Q 2 ^ip) = f(Al,A2); 



(15) 



A-independent pomeron flux and a parent nucleus dependent pomeron structure. In a 
model of this type ]3IJ the ratio Rai,A2 can be expressed using the nuclear structure 
functions F 2! a(%,Q 2 ) measured in inclusive electron nucleus scattering 



(16) 



2.4.3 Measurement of the A-dependence of the Fraction of Rapidity Gap Events 

One of the simplest measurements which could discriminate between the two pictures of pomeron 
formation proposed in |39|, f5)| and unresolved by the diffract ive ep scattering data, is the mea- 
surement of the A-dependence of the fraction of the number of rapidity gap events with respect 
to the total number of deep inelastic scattering events: 



R gap (Al,A2) 



N gap {Al)/N tot {Al) 
N gap (A2)/N tot (A2Y 



(17) 



In the model of |59] the pomeron is replaced by multiple soft colour exchanges between the 
quark- ant iquark pair into which the virtual photon has fluctuated and the target nucleus. In 
this model the ratio R gap (Al, A2) is expected to be 1. This is in contrast to the prediction 
of the colour singlet exchange model |4(| in which the ratio R gap (Al, A2) can reach a value of 
3 between nuclei with atomic numbers Al and A2 differing by 200. In order to illustrate the 
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Figure 11: The integrated ratio R gap (Al, A2) plotted as a function of xp. The error bars 
correspond to a luminosity of 1/A1 and 1/A2 pb -1 for each nucleus. 



statistical precision of such a measurement we show in Fig. |TT] the values of R gap (Al, A2) and 
their statistical errors. 

The statistical precision of R gap (Al, A2) is dominated by the statistical precision of N gap (Al) 



and N gap (A2) and was determined with help of the RAPGAP [38[ Monte Carlo by counting 



events generated within the kinematic domain defined by the following boundaries: Q > 2 
GeV 2 , \t\ < 0.5 GeV 2 and (3 > 0.05. The ratio shown in Fig. |Tl] was set to unity as we have 
not tried to model the nuclear dependence of the R gap (Al, A2). We observe that statistical 
precisions of ~ 2 % can be reached by collecting integrated luminosity of 1/A pb _1 for each 
nucleus. The systematic errors, as for the measurement of F^\f3,Q 2 ,X]p) are expected to 
be smaller than the statistical errors if two nuclei are stored simultaneously at HERA. This 
measurement, even if carried out for two light nuclei, can easily rule out one of the two models 
of pomeron formation. 



2.4.4 Elastic Vector Meson Production 

We present the accuracy expected for exclusive photoproduction of J/ij) mesons at Q 2 < 4 GeV 2 
(eA — > eAJ/ip) and exclusive production of p° mesons at Q 2 > 4 GeV 2 (eA eAp°). 

As discussed in detail later, elastic (or exclusive) production of vector mesons in the reaction 
ep epV, where V is a vector meson u>, 0, J/ip...), is thought to be of a diffractive nature. 
Recent calculations (cf. e.g. ET], S3, H41 ffH]) indicate that the cross section for these processes 



may depend on the gluon momentum distribution xG(x) probed at a value of x, the fraction 
of the nucleon's momentum carried by the gluon, x ~ (Q 2 + rriy)/W 2 , where my is the vector 
meson mass and W is the photon-nucleon centre-of-mass energy. 

We have determined the statistical accuracy with which the ratios 
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1 da A /dt\ t=0 



(18) 



A 2 da D /dt\ t=0 

can be determined for 4 He, 12 C, 32 S and 208 Pb. Fig. [l^ shows the results for J ftp production 
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Figure 12: Expected statistical accuracy (vertical error bars) on the ratio Ry for elastic 
photoproduction (Q 2 < 4 GeV 2 ) of J/if> mesons. In all cases a luminosity of 1 pb _1 /v4 for each 
nucleus was assumed. The horizontal bars indicate the size of the bins. 



We used DIPSI f46j, a Monte Carlo generator based on that describes the available 
ZEUS data on J /if) photoproduction. Events were generated in the range Q 2 < 4 GeV 2 and 
30 < W < 300 GeV; the J /if) was assumed to decay into e + e~ or fi + fi~ pairs. Events in which 
the decay leptons were outside the coverage of the barrel and rear tracking detectors of HI and 
ZEUS (approximately 34° < d < 164°) were rejected. A further reduction in the number of 
accepted events by a factor 0.2 was applied to account for efficiency and acceptance effects; this 
factor was taken to be independent of x. 

Only statistical uncertainties are shown. Possible sources of systematic uncertainties are 
the luminosity, the branching ratio, the global acceptance (including trigger and reconstruction 
efficiency, muon or electron identification etc.), the contamination from incoherent events, the 
feed- in from inelastic J /if} production (a la photon-gluon fusion) and the feed-in from if) 1 pro- 
duction. All the above contributions would largely cancel in a ratio for simultaneously stored 
nuclei, with the partial exception of the luminosity. In practice, for an integrated luminosity of 
1 pb _1 , the statistical uncertainty is expected to dominate. 

Fig. p~3| shows the expected statistical uncertainty for p° production at Q 2 > 4 GeV 2 . Here as 
well we used DIPSI, and the parameters were chosen so that the generator describes the ZEUS 



data ||7|. Events were generated in the range 4 < Q 2 < 100 GeV 2 and 30 < W < 300 GeV. 
The p° was assumed to decay into tt + tt~ pairs; pairs with masses between 0.3 and 1.4 GeV were 
considered. Events in which the decay pions were outside the coverage of the tracking detectors 
of HI and ZEUS (approximately 15° < •& < 164°) were rejected. A further reduction in the 
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Figure 13: Expected statistical accuracy (vertical error bars) on the ratio Ry for elastic 
production (Q 2 > 4 GeV 2 ) of p mesons. In all cases a luminosity of 1 pb~ 1 /A for each nucleus 
was assumed. The horizontal bars indicate the size of the bins. 



number of accepted events by a factor 0.6 was applied to account for efficiency and acceptance 
effects; this factor was independent of x. 

One can see from the figures that the accuracy of the measurements would be sufficient to 
discriminate between the colour transparency expectation of a ratio close to unity and the vector 
meson dominance expectation of the ratio Ry ~ A~ 2 ^ 3 ; cf. also the discussion in section pT6. 



2.5 Parton Propagation in Nuclei 



This has been studied in various fixed target experiments in the past |48], [49| . It is investigated 
by measuring the final state hadrons in the deep inelastic scattering. Since there will be many 
more than one hadron per event in experiments at HERA the statistical errors are not expected 
to be a limitation and high accuracy should be possible. This subject is considered in more 
detail in section [3]7| and in |H| . 



2.6 Other Physics 

In addition to the primary programme for nuclear beams in HERA outlined above there will 
be additional experiments which have not been explored in these studies. For example, there is 
physics interest in studying photoproduction from nuclear targets in the HERA energy range 
as well as sensitivity to physics beyond the standard model. Such sensitivity arises in e-nucleus 
collisions via coherent two photon interactions allowing production of positive C parity states 



which will be inhibited at LEP. This is discussed by Krawczyk and Levtchenko [51 
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3 Theoretical Overview 



3.1 Introduction 

In this overview we describe those aspects of the phenomenology of QCD which lead to the 
nonlinear effects referred to earlier. These effects are thought to be related to the mechanisms 
which will eventually limit the growth of the nucleon structure function as x decreases at finite 

Q 2 . 

3.2 Space-time Picture of DIS off Nuclei at Small x 
3.2.1 The Rest Frame 

In the rest frame of the target nucleus the life-time of a fluctuation is given by the formula 

r = ^— } (19) 

m N x Bj 

where (3 = Q 2 / (Q 2 + M 2 ) < 1. M is the mass of the qq system 

„ k? + m 2 

= i^iy m 

where z is the light-cone momentum fraction, k t the transverse momentum and m q the mass of 
the quark. Perturbative QCD studies show that the most probable configurations are those for 
which M 2 Q 2 . In the case of transversely polarised photons both configurations with small k t 
and highly asymmetric fractions z and configurations with comparable z and 1 — z contribute 
to the cross section. For the case of the longitudinal photons the asymmetric contribution is 
strongly suppressed. 

In the language of noncovariant diagrams this corresponds to the virtual photon fluctuating 
into a quark-antiquark pair at a longitudinal distance l c = from the nucleus which far 
exceeds the nuclear radius. The distance l c is referred to as the "coherence length". The pair 
propagates essentially without transverse expansion until it reaches the target. QCD evolution 
leads to a logarithmic decrease of (3 with increasing Q 2 . At HERA coherence lengths of up to 
1000 fm are possible, so that the interaction of the qq pair with nuclear matter can be studied 
in detail - notably its transparency to small size pairs - colour transparency. 

At HERA new features of colour transparency should emerge: the incident small size qq 
pair resolves small x gluon fields with virtualities ~ Q 2 . If the transverse size of the qq pair is 
r t = b q — bq , the cross section for interaction with a nucleon is [52 



vr 2 

o- qq -, N (E inc ) = —r 2 a s (Q 2 )xg N (x,Q 2 ), (21) 

where Q 2 ~ X(x ~ 10 -3 ) ~ 9, x = ^n^E — ■ Since the gluon density increases rapidly with 
decreasing x, even small size pairs may interact strongly, leading to some sort of perturbative 
colour opacity - the interaction of a small object with a large object with a cross section 



comparable to the geometric size of the larger object (Fig. |14]). 

Unitarity considerations for the scattering of a small size system [13] - i.e. the requirement 
that ai ne i(qq,target) < ^R 2 arget - indicate that nonlinear effects (i.e. effects not accounted for 
by the standard evolution equations) should become significant at much larger values of x in 
eA scattering than in ep scattering. 
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Figure 14: Colour-dipole cross section, UqqN^x, b) of Eq. (0), function of the transverse 
size of the qq pair for various values of x and for the GRV94 parametrization of the nucleon's 
gluon density. 



3.2.2 The Breit Frame 

In the Breit frame, small x partons in a nucleon are localized within a longitudinal distance 
~ 1/xpN, while the distance between two nucleons is ~ r NN m N /p N (r NN is the distance 
between nucleons in the rest frame and p^ is the nucleon momentum). Therefore partons with 
x < 1/(2771^^), where ta ~ tqA 1 ^ 3 fm is the nuclear radius and r$ = 1.1 fm cannot be localized 
to better than the whole nuclear longitudinal size. Hence low x partons emitted by different 
nucleons in a nucleus can overlap spatially and fuse, provided the density is high enough, 
leading to shadowing of the partonic distributions in bound nucleons with respect to the free 
nucleon ones and to nonlinear effects already at values of x ~ 10~ 4 4- 10 -3 . For example, in the 



simplest model of nonlinear effects corresponding to the fan diagrams of Fig. [15], the additional 
contribution Sg^ix, Q 2 ) to gA(%,Q 2 ) due to the nonlinear term in the equation for the Q 2 ,x 
evolution of the gluon density is fl9fl : 



2 d 6xg A (x,Q 2 ) 81 A 1 / 3 2 2 r 1 du , - 

Q W j = -T6Q^ a ' iQ } / - U9nM } • (22) 



u 



The analogous equation for the gluon density in the nucleon has a much smaller coefficient - 
approximately by a factor r^/rf^A 1 / 3 , where ~ 0.8 fm is the nucleon radius. Once again 
one can see then that the x-range where nonlinear effects become significant differs for a heavy 
nucleus and for a nucleon by more than two orders of magnitude, assuming xg^(x,Q 2 ) oc x n 
with n ~ —0.2. 

Thus electron-nucleus collisions at HERA can be seen as efficient amplifiers of nonlinear 
QCD effects. 



3.3 Theoretical Framework for Small x Phenomena in eA Collisions 
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Figure 15: Typical fan diagrams leading to nonlinear evolution of g^ix^Q 7 



3.3.1 Perturbative and Nonperturbative Shadowing 



At small x the DIS cross section per nucleon in a nucleus is smaller for a bound nucleon than for 
a free one, the so called shadowing phenomenon. Shadowing is determined by a combination of 
non-perturbative and perturbative effects. In the DGLAP evolution equations one can express 
shadowing at large Q 2 through the shadowing at the normalization point Qq. This type of 
shadowing is connected to the soft physics. It can be visualized e.g. in the aligned jet model of 
Bjorken 0, extended to account for QCD evolution effects |5B| . A virtual photon converts to a 
qq pair with small transverse momenta (large transverse size) which interacts with the nucleus 
with a hadronic cross section, leading to shadowing. The effective small phase volume of these 
configurations (oc leads to Bjorken scaling and it is due to colour transparency p3 |. 

At large Q 2 , these qq pairs evolve into systems with gluons, leading to a shift of shadowing 
to smaller x, which is equivalent to the standard Q 2 evolution of parton distributions. These qq 
pairs, which interact with the target nonperturbatively, seem to be responsible for most of the 
shadowing at intermediate Q 2 and x ~ 1CT 2 which has been studied at fixed target energies. 
This mechanism of shadowing is effective for <jt only since for <tl the aligned jet contribution 
is strongly suppressed. For (as well as for the production of heavy quarks) one is more 
sensitive to the shadowing due to the interaction of small size qq pairs with the nuclear gluon 
field which can be shadowed. 

At smaller x the situation may change rather dramatically because, as the recent HERA data 
indicate, already for Q 2 ~1.5 GeV 2 at x ~ 1CT 4 perturbative contributions to F 2p (x, Q 2 ) appear 
to become important, leading to a rapid increase of the structure functions with decreasing x. 
Hence contributions of various perturbative mechanisms which may generate shadowing for 
configurations of a size smaller than the hadronic size may become important. Perturbative 
QCD may be applicable to those small size pairs. Typical contributions involve diagrams of 
the eikonal type, various enhanced diagrams, etc. (Fig. p~5| . p~6| ) - 



3.3.2 Shadowing and Diffraction 

In practically all models it is assumed that nuclei are built of nucleons. So the condition that 
the matrix element < A\T[J ll (y)J u (0)}\A > involves only nucleonic initial and final states is 
implemented^. Under these natural assumptions one is essentially not sensitive to any details 

1 The condition that nuclei are built of nucleons is not so obvious in the fast frame picture. However it is 



implemented in most of the models |19 
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Figure 16: Examples of typical perturbative QCD diagrams contributing to nuclear shadowing. 



of the nuclear structure, such as short-range correlations etc. 

In the case of scattering off the deuteron and light nuclei the same diagrams contribute to 
the cross section for diffraction in ep scattering and the cross section for shadowing - hence 
similar nonlinear phenomena like those described by eq. (|22]) are involved in each case. For 
example for the deuteron M: 



0~shad 



a tot {eD) - 2a tot (eN) 



dt \t=0 



where R 



(l-A 2 ) 



A 



a(eN) 
Re A/ ImA « 7vdlnA 



crtot(ep) 



8nR 2 



■R, 



(23) 



D 



.-, for the amplitude A of 7*p scattering and Rjj 
is the deuteron radius. For small x, A may be as large as 0.5, leading to R ~ 0.5 especially 
for the case of the longitudinal cross section. So already for light nuclei the study of the 
total cross sections of scattering from nuclei would allow to establish a fundamental connection 
between the two seemingly unrelated phenomena of diffraction at small t in ep scattering and 
nuclear shadowing. With the increase of A more complicated nonlinear interactions with several 
nucleons become important, see e.g. Fig. |TB|b. 

Nuclear shadowing for the total cross sections has a simple physical meaning - it corresponds 
to a reduction of cross section due to screening of one nucleon by another (as well as by several 
nucleons for A > 2). If one treats the deuteron as a two nucleon system it is possible to apply 
the Abramovskii, Gribov, Kancheli (AGK) cutting rules |55| to elucidate the connection 
between nuclear shadowing, diffraction and fluctuations of multiplicity. One observes 
that the simultaneous interaction of the 7* with the two nucleons of the deuteron modifies not 
only the total cross section but also the composition of the produced final states. It increases 
the cross section for diffractive scattering off the deuteron due to diffractive scattering off both 
nucleons by Sadiff = <J s h a d- At the same time the probability to interact inelastically with one 
nucleon only is reduced since the second nucleon screens the first one: Sa S i ng i e = —Aa s h a d- In 
addition, a new process emerges in the case of the deuteron which was absent in the case of 
the free nucleon - simultaneous inelastic interaction with both nucleons which leads to a factor 
of two larger multiplicity densities for rapidities away from the current fragmentation region: 
^double = ^shad- Altogether these contributions constitute —o~ s had, the amount by which the 
total cross section is reduced f\ 



2 For simplicity we give here relations for the case of purely imaginary 7*iV amplitude -5^4 = 0. 
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To summarize, there is a deep connection between the phenomena of diffraction observed at 
HERA in ep scattering and nuclear shadowing as well as the A-dependence of diffraction and 
the distribution of the multiplicities in DIS. 

It follows from the above discussion that it is possible to get information about the dy- 
namics of nuclear shadowing and hence about nonlinear effects by studying several key DIS 
phenomena such as: nuclear shadowing for inclusive cross sections F^, pA chaTm ^ the cross 
section for nuclear diffraction; the multiplicity distribution for particle production in the central 
rapidity range; diffractive production of vector mesons. The advantage of the latter process is 
that one gets a rather direct access to the interaction of a small colour dipole with matter. It 
is in a sense an exclusive analogue of o~l which is easier to measure. 



3.4 The ^-Dependence of Parton Distributions at Small x 

As discussed above, the nucleus serves as an amplifier for nonlinear phenomena expected in 
QCD at small x. The simplest example of such effects is given by equation (|22"D where the 
nonlinear term is proportional to the square of the nucleon gluon density. If shadowing were 
absent the parton densities per unit transverse area would be enhanced by a factor A 1 ^ as 
compared to the free nucleon case. Hence even just an upper limit on the parton densities 
based on unitarity - that the cross section for the inelastic interaction of a small dipole with 
a nucleus may not exceed Oi ne \ = 7rR\ - leads to the expectation of nonlinear phenomena - 



shadowing of an observable magnitude - already at a; ~ 10 4-10 [|4|. 

Hence, from detailed studies of the A-dependence of the parton densities it would be possible 
both to check the dominance of the two-nucleon screening mechanism for x ~ 10 -2 J56|, [5?J and 
to extract information about the coherent interaction of the virtual photon with three (four) 
nucleons at x < 10~ 3 . 

For x > 10~ 2 for any nucleus and for all x for light nuclei, the main contribution to shadowing 
is given by the interaction with two nucleons of the target. Hence in this regime there is a 
relatively simple connection with the diffraction of a virtual photon off a proton - which is the 
simplest nonlinear effect in the perturbative domain in QCD. For smaller x and heavy nuclei, 
when essential longitudinal distances become comparable and ultimately exceed the diameter 
of the nucleus, several nucleons at the same impact parameter contribute to the screening. 
It is worth emphasizing that these multi-vacuum exchange processes cannot be singled out 
unambiguously using a nucleon target. The relevant QCD diagrams for the total cross section 
of 7* A interaction are rather similar to higher-order nonlinear diagrams for the proton target - 
except that in the nuclear case one has to impose the condition that couplings to the individual 



nucleons are colour singlets, see e.g. Fig. 15,16 



In a sense, the studies of nuclear shadowing at small x and large Q 2 can be considered as a 
simpler model of nonlinear effects which occur in the case of a nucleon target. In the latter case 
it is not easy to relate the coupling of say two vacuum exchanges (or a ladder with 4 gluons 
in the i-channel) with a nucleon to the coupling of one vacuum exchange with a nucleon. In 
fact the region of 10~ 3 > x > 10~ 4 may be optimal in this respect since nonlinearities for 
the nucleon case are still small though nonlinearities for the nuclear case are already quite 
substantial. It is worth emphasizing that experience of the studies of the total hadron-nucleus 
scattering indicates that interaction with bound nucleons for the total cross sections can well be 
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Figure 17: A- dependence of nuclear shadowing and probability of rapidity gap events in the 
colour screening model of shadowing; dot-dashed curve assumes A-independent probability of 
rapidity gap events. 



approximated by the interactions with free nucleons (for a recent analysis see Therefore 
nuclear structure effects do not obscure the interpretation of nuclear shadowing effects. 

Using current information from HERA on diffractive production in ep scattering it is 
straightforward to estimate the amount of nuclear shadowing at small x taking into account 
interactions with 3 or more nucleons using the eikonal approximation with an effective cross 
section determined from diffractive data, see eq.(^) below. The result of the calculation [[KJ 
is shown in Fig. [17] for Re/Im = 0; for A > 12 it weakly depends on the value of Re/Im. 
Since the data on diffraction indicate that the fraction of diffractive events in DIS weakly de- 
pends on x, Q 2 these considerations show that significant shadowing effects should be present 
for FA(x, Q 2 ) in the wide small x range of HERA. Note that the shadowing effect in DIS is 
expected to be much smaller than for the case of real photon scattering since the effective cross 
section for interaction of the hadron component of quasi-real photon at HERA is a factor of 
~ 3 larger than for a highly virtual photon (we use here the HERA data on diffraction for real 
photons floT 



Since the interaction of the octet colour dipole gg is a factor of 9/4 stronger than for the 
qq dipole, nonlinear effects are expected to be more important for gluons. So gluon shadowing 
would provide even more direct access to nonlinear phenomena. Note that in this case there is 
no simple relation of shadowing with diffraction in 7* +p DIS, so any information about gluon 
shadowing would be complementary to the information from ep DIS. There are very few data 
on the gluon distribution in nuclei. Among them, the enhancement of the gluon distribution 
at x ~ 0.1 indicated by the inelastic J/ip production data [25|. Also the analysis |j0] of 
the scaling violation for the ratio F^ n jF^ |TT| under the assumption that higher twist effects 
are not important in the Q 2 ,x range of the data allows to extract information about the A- 
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dependence of gluon distributions, indicating some nuclear shadowing for G a for x < 0.01 and 
an enhancement at x ~ 0.1, see Figure in |]6O|0- Theoretical expectations for gluon shadowing 
discussed in the literature are quite different - from a larger effect than for FA |63[ , to an effect 



comparable to that of quarks |)2| [61], |64], Q to substantially smaller shadowing |p6| ; see also 
contributed papers to these proceedings. 

Comparison of different determinations of shadowing of gluons and measurements of the 
scaling violation for the F^/F® ratios will allow to determine the range of applicability of the 
DGLAP evolution equations and hence provide unique clues to the role of nonlinear effects. 

It is worth emphasizing also that knowledge of parton distributions in bound nucleons at 
these values of x will be crucial also for studies in heavy- ion physics at the LHC and RHIC. 



3.4.1 BFKL Pomeron 

One can envision several strategies for the study of the BFKL Pomeron in DIS. The main 
requirement is to enhance the contribution of scattering of small transverse size configurations 
in the ladder. It is natural to expect that screening for these configurations would be minimal. 
Hence for heavy nuclei the contribution of the BFKL Pomeron can be enhanced. 

1. A procedure can be envisioned to study the A-dependence of F^(x,Q 2 ) at small x to 
extract the term in the structure functions oc A and then study its x dependence. Based 
on the above argument A.Mueller has suggested [|67| that the x dependence of this term 



(linear in A) would be closer to BFKL type behaviour. 

2. One promising direction to look for the BFKL pomeron is the Mueller-Navalet process of 
producing two high p t jets with large a rapidity difference p)8[] to suppress the contribution 
from small transverse momenta (large transverse distances) in the ladder. In the case of 
a nuclear target large distance contributions would be screened out to a large extent. 

3. Another possibility is the production of p mesons at large \t\ in inelastic diffraction. To 
enhance the contribution of the BFKL Pomeron it is desirable to increase the contribution 
of small configurations in the p meson, i.e. quark-antiquark pairs with small transverse 
separation. Large size configurations can be filtered out by exploiting the fact that they 
are absorbed on the nucleus surface. Once again extracting the term in the cross section 
oc A would allow to enhance the contribution of the BFKL Pomeron. 



3.5 Diffraction off Nuclei 
3.5.1 Introduction 

Diffraction studies have been defined as one of the primary goals of nuclear beams in HERA. 
Such processes can be interpreted using two complementary languages depending on whether 
the rest frame of the nucleon or the Breit frame are used: 



Scattering of electrons on colourless components of the proton 36, p%. Such scattering 



may be identified, for the very low x events dominated by diffraction, with the interaction with 



3 The shadowing for gluons should be accompanied by a significant enhancement at larger x since the total 
momentum fraction carried by gluons in nuclei is not suppressed and is probably slightly enhanced fl3l| . 
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the vacuum t-channel exchange which is often referred to as the Pomeron, IP. This object 
is not necessarily the same as the Pomeron of the Gribov-Regge high-energy soft interactions 
(see report of the diffractive group). Deep inelastic electron scattering leading to the presence 
of a rapidity gap can thus be considered as probing the internal parton structure of the IP 
originating from the proton. 

One of the questions of primordial importance which may be addressed within the future 
electron-nucleus scattering program at HERA is then "how universal is the internal structure 
of the Pomeron?" or, more precisely: "Is the internal structure of the Pomeron originating 
from various hadronic sources (protons, neutrons, nuclei) the same?". We shall show below 
how nuclei may help in answering these questions. 

• The diffractive interaction of different hadronic components of the virtual photon with 
the target via vacuum exchange. Diffraction predominantly selects the 7* components which 
interact with sufficiently large cross sections such as large transverse size qq, qqg colour dipoles. 
Therefore the study of diffraction plays a very important role in determining the relative impor- 
tance of small and large size configurations and addressing the question whether small white 
objects interact weakly or not. Indeed if the interaction with a target becomes sufficiently 
strong at small impact parameters the cross section for diffraction (which includes both elastic 
scattering and inelastic diffractive dissociation) would reach the black body limit of 50% of the 
total cross section. 



3.5.2 Theoretical Expectations 

Diffraction off a nucleon (including dissociation of the nucleon) constitutes about 15-20% of 
the deep inelastic events. Therefore the interaction is definitely far from being close to the 
scattering off a black body. Even this number came a surprise in view of the large Q 2 value 
involved. Using the generalized optical theorem as formulated by Miettinen and Pumplin, one 
can estimate the effective total cross section for the interaction of the hadronic components of 
the 7* as 

a eff = 16vr -. — « 12 + 15mb. (24) 

This cross section is significantly smaller than the pN cross section which at the HERA energies 
can be estimated to be close to 40 mb using the vector dominance model and the Landshoff- 



Donnachie fit [69]: 



= °&(8 )(±) , (25) 
Vs / 

where n ~ 0.08, so = 200 GeV 2 , crfj^(s ) = 25 mb. However it is sufficiently large to result in a 
substantial cross section of diffraction for small x - it can reach 30-40% for large A (Fig. |i~7|) ||40|| . 
For large A the coherent diffraction dominates when the incoming wave is sufficiently absorbed 
at small impact parameters which, by virtue of Babinet's principle, corresponds to scattering 
beyond the nucleus. In such processes the nucleus remains intact and the average momentum 
transfer is very small ((£) oc A~ 2 ^ 3 ). 

One expects that hadronic configurations interacting with different strength contribute to 
diffraction (cf. Fig. |1^). The parameter a e ff characterizes just the average strength of this 
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interaction, while the distribution over the strengths is expected to be quite broad. The study 
of diffraction off nuclei allows to separate contributions to diffraction of large and small size 
configurations due to the filtering phenomenon: with the increase of A the relative contri- 
bution of more weakly interacting (smaller size) configurations should increase since they are 
less shadowed, leading to a relative enhancement of the colour transparent subprocesses. 

Examples of promising processes are: 

• Diffractive production of charm. The A-dependence of this process would be interesting 
already at low Q 2 since the essential transverse distances are, naively, of the order of 
f/m c , where m c is the charm quark mass. Since the cross section for the interaction 
of a colour dipole of such size is small for x ~ I0~ 2 , the cross section for diffractive 
charm production at these values of x is small and practically not shadowed, leading to 
a cross section oc A 4 / 3 . At the same time the cc — N cross section increases rapidly with 
decreasing x (increase of energy) for fixed Q 2 . Therefore at HERA energies diffractive 
charm production in ep collisions may become a substantial part of the total diffractive 
cross section. At this point one expects the emergence of shadowing in diffractive charm 
production in eA collisions, leading to slowing down of the A-dependence of diffractive 
charm production as compared with the A dependence at x ~ I0~ 2 . 

One can go one step further and study the A-dependence of p t distributions for diffractive 
charm production. Smaller size components will be less absorbed and so their relative 
contribution may increase with A. 

• Diffractive production of two high p t jets. 

Selection of large p t jets enhances the contribution of diffraction of small size configura- 
tions. Hence, one expects broader p t distributions in the case of nuclear targets (smaller 
jet alignment) with nontrivial dependences on W and Q 2 . For example, if we fix the p t 
of the jets, the A-dependence of dijet production should become weaker with increasing 
energy reflecting the increase of the absorption (which can be studied this way). If on 
the other hand we fix W and consider the A-dependence as a function of pt, a stronger 
A-dependence is expected. Effectively, this would be another way to approach colour 
transparency via filtering out of the soft components. 

To summarize, a study of inclusive diffraction will give better insights into the structure 
of the Pomeron. The interplay of soft and hard contributions will lead to a breakdown of 
factorization for the structure function of the Pomeron. Stated differently, the check of the 
degree of "universality" of the Pomeron - i.e. whether the "nuclear" Pomeron is different from 
the Pomeron observed in ep diffraction - will provide a very sensitive test of QCD dynamics. 

An important aspect of the diffractive studies is the colour transparency phenomenon. In 
view of its special interest we will discuss this separately below. 



3.5.3 Multiplicity Fluctuations 

As we explained in section |3.3.2| diffraction off nuclei and nuclear shadowing are closely related 
to the simultaneous inelastic interactions of the virtual photon with several nucleons. Such in- 
teractions produce events with large multiplicity densities in the central rapidity range, leading 
to a much broader distribution over multiplicities in eA collisions than in the ep case. Study 



26 



The ratio of charged particle multiplicities distributions 
for 80 < W < 115 GeV and < eta < 2 for two nuclei 



[dNch(A1 )/dn)/(dNch(A2)/dn) 

statistical errors for the luminosities 1 /A pb— 1 

-<J>- - color singlet model A1 =200, A2=1 
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Figure 18: A-dependence of distribution over multiplicity calculated in the colour singlet model 
for a luminosity of 1 pb -1 . 



of this effect will provide information complementary to that obtained from nuclear shadowing 
about the structure of the vacuum exchange at small x. Interesting phenomena to look for may 
be: 

1. Local fluctuations of multiplicity in the central rapidity region, e.g. the observation of 
a broader distribution of the number of particles per unit rapidity, n(Ar]), than in ep 



scattering IptOfl , see Fig. 18 



2. Long range rapidity fluctuations - i.e. positive correlation of the increase of multiplicity 
in one rapidity interval with the increase of multiplicity several units away. 

3. Correlation of the central multiplicity with the multiplicity of neutrons in the neutron 
detector (most effective for heavy nuclei). 



3.6 Colour Transparency Phenomena 

An important property of QCD is that small objects are expected to interact with hadrons 



with small cross section ||70|| . This implies that in the processes dominated by the scatter- 
ing/production of hadrons in "point-like" (small size) configurations (PLC) when embedded in 
the nuclei, the projectile or the outgoing hadron essentially does not interact with the nuclear 
environment J7T], [7Hfl . In the limit of colour transparency one expects for an incoherent cross 
section a linear dependence on A, for example 

da{e + A^e+p+{A- 1)*) _ da(e + p -> e + p) 
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Figure 19: Average transverse size of the qq components effective in A^n->vn for p— and 
J/^-meson production and ctl- The probed Q 2 scale is inversely proportional to b 2 . 



while for coherent processes at t = one expects 



da(Y + A -> X + A) 
dt 



A 



2 da(Y + N 



X + N) 



dt 



(27) 



No decisive experimental tests of this property of QCD were performed so far since in most of the 
current experiments the energies were not sufficiently high to prevent expansion of the produced 
small system. The high-energy E665 experiment [7^] at FNAL has found some evidence for 
colour transparency in the p meson production off nuclei. However, the data have low statistics, 
cover a small x, Q 2 range and cannot reliably separate events without hadron production. 

A quantitative formulation of colour transparency for high-energy processes can be based 



on eq.(21). For the case of nuclear targets it implies that for a small enough colour dipole, 
the cross section of its interaction with nuclei is proportional to A up to the gluon shadowing 
factor. As a result the colour transparency prediction for 2 jet and vector meson diffractive 
production is J73], f|2| Q 
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(28) 



Gluon shadowing constitutes a rather small effect for x ~ 10~ 2 (see earlier discussion). For 
smaller x it increases but it is in any case much smaller than the screening effect expected 
in the case of lack of colour transparency if the produced system interacts with cross section 
comparable to <j p n ~ 30-40 mb. For such values of o one expects the cross section to behave 
as oc A 4 / 3 for t — which would be possible to test using diffractive production by quasi-real 
photons. 

Coherent diffractive p, J -meson production 



4 In writing eq. (|2g) we neglect the difference of Q 2 scales for different processes which is reflected in a different 
dependence of the essential transverse size of the qq state on the process (see Fig. 19). For a discussion of the 
appropriate scale for dijet production see ]7q|. 
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The most straightforward test of colour transparency can be made using coherent production 
of p or J/^-mesons at small t using nuclei with A > 12. The p t resolution of the current 
detectors is good enough to single out the diffractive peak which is concentrated at p t < 0.1 
GeV. In the higher x end of the range which could be studied at HERA for vector meson 
production, x ~ 10~ 2 , one expects at large Q 2 nearly complete colour transparency since gluon 
shadowing effects are rather small and decrease rapidly with increase of Q 2 , while the transverse 
separation, b, between q and q is of the order of 0.4 fm for Q 2 ~ 10 GeV 2 and further decreases 
with increase of Q 2 (Fig. [19| ). Study of coherent J/ip meson production would allow to 
probe colour transparency for propagation of even smaller dipoles since {b c5 (Q 2 = 0)) ~ 0.2 fm. 

On the other hand as discussed earlier at the smallest values of x of the HERA range, 
screening effects should start to play a role even at large Q 2 so a gradual disappearance of 
colour transparency is expected - the emergence of colour opacity. Noticeable screening is 
expected already on the basis of unitarity constraints. Qualitatively one may expect that the 
rise of the cross section for vector meson production with increasing energy at fixed Q 2 will slow 
down at significantly lower energies than for the case of the 7*+p reaction. Currently theoretical 
calculations of vector meson production by transversely polarised photons are difficult because 
the nonperturbative large distance contribution is not as strongly suppressed in this case as in 
the longitudinal case. If contribution of pairs with large transverse size is indeed important for 
<tt, it would be filtered out with increasing A leading to larger values of <j l /<jt for large A. 

Let us enumerate several other effects of Colour Transparency (CT) in diffractive production. 



1. Production of excited vector meson states p', <ft'. 

In the CT limit, QCD predicts a universal A-dependence of the yields of the lowest mass 
and excited states (this includes the effect of gluon shadowing in eq.(|28|)). This is highly 
nontrivial since the sizes of the excited states are much larger, so one might expect larger 
absorption. On the other hand, for lower Q 2 average transverse distances, b, are not 
small. At these distances the wave functions of ground and excited states differ. So in 
this Q 2 range the relative yields of various mesons may depend on A. 

2. Production of high p t dijets. 

The uncertainty relation indicates that coherent production of dijets with large p t , car- 
rying all the momentum of the diffractively produced system is dominated by distances 
r t oc Hence filtering out of soft jets is expected, leading to a broadening of the p t and 
thrust distributions. At the same time the study of the A-dependence of low p t jets would 
allow to address the question of colour opacity. Another feature to look for would be the 
distribution over the electron-two jet plane angle as suggested in J75fl. 

3. Coherent diffractive production at —t > 0.1 GeV 2 for A = 2,4. 

An important question here is the possibility to observe the "disappearance" of colour 
transparency for p-meson production and the emergence of "colour opacity" - due to 
nonlinear screening effects at x ~ 10 -4 . Manifestations of CT would be the increase of 
the differential cross section ^ below the diffractive minimum (|i m j n ( 4 i7e)| ~ 0.2 GeV 2 
and suppression of the cross section in the region of the secondary maximum. A gradual 
disappearance of CT in this region with increasing energy would appear as a very fast 
increase with energy of the secondary maximum of the t distribution. Remarkably, in this 
region the cross section for the process is proportional to [Gn(x, Q 2 )} 4 , where G^ is the 
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gluon density in the nucleon [77, 78}. The present beam optics would allow measurements 
of quasielastic processes with 4 He in the region of the secondary maximum (|t( 4 He)| pa 
0.4 GeV 2 ). For a luminosity of 10 pb _1 it would be possible to measure the p-meson 
production cross section up to Q 2 ~ 10 GeV 2 . 

4. A-dependence of rapidity gaps between jets in photoproduction. 

Recently photoproduction events which have two or more jets have been observed in the 
range 135 < W 1P < 280 GeV with the ZEUS detector at HERA [[7|. A fraction of 
the events has little hadronic activity between the jets. The fraction of these events, 
/(At/), reaches a constant value of about 0.1 for large pseudorapidity intervals At/ > 
3. The observed number of events with a gap is larger than that expected on the basis of 
multiplicity fluctuations assuming the exchange of a colour singlet. This value is rather 
close to estimates in perturbative QCD 



80, 81 



to interactions of spectator partons 
and CDF 



S3 



81 



2) neglecting absorptive effects due 
It is much larger than the values reported by DO 
Small effects of absorption are by no means trivial in view of the 
large interaction cross section for many components of the hadronic wave function of the 
real photon. They may indicate that colour transparency is at work here as the ZEUS 
trigger may select point-like configurations in the photon wave function [85]. To check 



this idea it would be natural to study the A-dependence of rapidity gap survival. It is 



demonstrated in |85j that this probability strongly depends on the effective cross section 
of the interaction of the photon with the quark-gluon configurations involved in producing 
rapidity gap events. One would be sensitive to cross sections as small as ~ 5 mb. 



3.7 Parton Propagation in Nuclear Matter 
3.7.1 Introduction 

Measurements of final state hadrons allow to investigate the effects of partons propagating 
through nuclear matter. In the present section we discuss some of these possibilities. The 
discussion is restricted to incoherent phenomena (coherent nuclear interactions were discussed 
in previous sections). 

There are essentially two types of measurements which can be useful: (a) energy loss of high- 
energy particles and (b) increase of their transverse momentum, both studied as a function of 
the nuclear number A and/or number of nucleons emitted from the target nucleus. They are 
sensitive to different aspects of the interactions. Energy loss reflects the properties of inelastic 
collisions: the value of the inelastic parton-nucleon cross section and of the inelasticity. The 
increase of transverse momentum and emission of the nucleons from the target can be induced 
by elastic as well as by inelastic collisions and thus can provide information on both. 

At this point we emphasize the importance of the measurement of the distribution of the 
nucleons (protons and neutrons) emitted from the target nucleus during or after the interaction. 
Such measurements give direct access to the number of secondary interactions inside the target 



36 , as was already realized (and used) in numerous emulsion experiments where protons with 



momenta 250 < < 700 MeV/c were measured [57]. Related information can be inferred 



from the measurement of the production of soft neutrons (E n < 10 MeV) which was studied 
recently by the E665 collaboration |88|. Measurements of the emitted nucleons and of their 



energy spectrum should thus be considered a high priority. They allow to improve greatly the 
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analysis of the data in three respects: (a) by studies of the distribution of nucleons themselves 
one obtains information of the strength of the secondary interactions and - more essentially- on 
the fluctuations which are expected to be large and would otherwise hamper the interpretation 
of the data; (b) By studying the interactions as a function of the number of the emitted nucleons 
one can obtain not only a larger lever arm in terms of the number of secondary interactions 
but also, at a fixed nuclear number, a clean sample, free of possible biases related to the use of 
different targets; (c) a really exciting possibility is to develop a trigger for events with a large 
number of emitted nucleons (or highly reduced charge of the nuclear remnant). This would 
allow to study in detail the rare events with particularly strong secondary interactions. Further 
quantitative studies (both theoretical and experimental) of this problem are necessary to clarify 
the relation between the emitted nucleons and number of secondary interactions. Some work 
on these lines was already presented during this workshop |89|, |J |9(J. The detection of such 
nucleons at HERA will be simpler than in fixed target experiments since they will be boosted 
by the motion of the target nucleus ||. 

As we discussed previously there should be a substantial difference between the nuclear 
effects observed in the region of "finite" x (x greater than, say, 0.05) and the region of very 
small x (x smaller than, say, 0.001). The reason is the different nature of the photon-nucleon 
interactions in these two regions caused by the difference in life times of the relevant photon 
fluctuations into a qq pair. 

One sees from eq.( [19]) that at finite x the life time of a fluctuation is rather short (smaller 
than 2 fm). This has two consequences: (a) the interaction of the photon must take place inside 
one of the nucleons of the target (nuclear coherence suppressed) and (b) the high-energy part 
of the interacting photon fluctuation can be well approximated by a simple "bare" quark. This 
can be seen as follows. The time necessary to produce a high-energy gluon from a quark is 
given by 

2-Ev, 

r 9 = -f (29) 

Ht 

where E g is the gluon energy and q t its transverse momentum (with respect to the quark). 
From the obvious condition r g < r we then obtain 

E ' < 2^' < 3 °> 
We conclude that at finite x there is simply no time to produce the energetic gluon cloud. 

For DIS nuclear collisions this picture implies that after the first interaction of the virtual 
photon the nucleus is penetrated by one bare quark (following approximately the direction 
of the virtual photon). Alternatively one can consider the process of production of two high 
Pt jets in photon- gluon fusion. In this case one studies propagation through the nucleus of a 
colour octet state. Thus the observed nuclear effects measure interactions of a bare parton (or 
a system of bare partons) in nuclear matter. 

The situation is rather different in the region of small x. In this case there is enough time 
for the qq system coupled to the virtual photon to "dress" itself into a cloud of energetic gluons 
and qq pairs (in the limit xbj — » the condition ( p0|) is not restrictive). Consequently, the 
system traversing the nucleus is a complex multiparton system resembling in some respects 
an "ordinary" hadron. The observed nuclear effects measure interactions of this multiparton 
system {dressed quark or qq dipole) in nuclear matter. It should thus not be surprising that 
the expectations for the region of small x are rather different from those at finite x. 
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3.7.2 Current Experimental Situation 



Two major manifestations of the interaction of a parton propagating through the nuclear 
medium which were studied experimentally so far are the parton energy loss and the broadening 
of its transverse momentum distribution. 

The measurements of the leading hadron spectrum in deep inelastic lepton-nucleus scattering 
have the biggest sensitivity to the energy loss. Measurements at incident energies below 50- 
100 GeV find a depletion of the leading particle spectrum which may be interpreted as due 
to the energy loss. At higher energies the effect nearly disappears indicating that energy 



losses for partons propagating through nuclear matter are definitely smaller than AE/dz < 1 
GeV/fm. This is in agreement with the observed low multiplicity of low energy neutrons in 
fiPb interactions at high energies f88f , which is consistent with knockout of one nucleon from 



lead [Q. This makes it impossible to address directly the question of energy losses at high 
energies. 

The phenomenological situation with transverse momentum broadening is somewhat con- 
fusing at the moment. The //-pair production experiments |9l[] which measure the broadening 
of the incident quark find rather small broadening: Apf ~ 0.1 GeV 2 for the distance, L ~ 5 
fm. At the same time the transverse momentum broadening for the jets produced by outgo- 



ing partons in photon-nucleus interactions seems to be much larger [92]. Theoretically this 



difference is not understood [93 



3.7.3 Perturbative QCD Expectations for Finite x 

All estimates of nuclear effects in this region of x accept that the interaction of a bare quark 
in nuclear matter is dominated by colour-exchange processes which lead to break-up of the 
"wounded" nucleon in the target but do not slow down significantly the energetic quark. This 
is based on the argument that the quark in question can only emit gluons satisfying the condition 
T g < Al where Al is the distance between the subsequent collisions. From eq.(p9|) we deduce 
that the energy loss in one secondary collision is limited by 

AE < 2g 2 AZ, (31) 

i.e., by a value which is independent of the energy of the quark. Consequently, a high energy 
quark can lose only a small fraction of its energy. A precise evaluation of the energy loss in a 
collision with one target nucleon is not possible at the present stage of the theory. 

However recently a significant progress was obtained in the analysis of the propagation 
of a virtual parton through the nuclear medium [[|. It was demonstrated that the Landau- 
Pomeranchuk-Migdal effect in QCD is qualitatively different from the case of QED. It was 
argued that for sufficiently large distances, L, traversed by a parton the process is dominated 
by perturbative QCD though the momentum transfers in the individual collisions are small. 

A simple relation was found between the p t broadening and the energy loss 

dE a s N c 2 
-^ = — Ap,(L), (32) 

which corresponds to substantially smaller energy losses than those implied by the inequality 
in eq. (|3l|) . 
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Probably the most striking prediction is the quadratic dependence of the energy loss on the 
traversed distance, L, as compared to the nonperturbative models where it is approximately 
proportional to L. Numerically the authors find for a quark 

- Arart 1nTfa) 2 ' (33 > 

neglecting the x dependence of the nucleon gluon density (and a factor of ~ 2 larger energy 
loss for gluons). If this effect is included the L-dependence is even steeper. The numerical 
coefficient is estimated here from the information on the nucleon gluon densities and consistent 
with eq.( j32|) if one uses experimental data on p t broadening of the /z-pair spectrum |9l[ 



Alternatively, if one uses the p t information from 7 A data a much larger energy loss is predicted. 
However even in this case the expected energy loss is too small to be observed directly in DIS 
at collider energies. 

The energy loss occurs (for realistic nuclei) via the emission of one or two gluons with 
energies ~ AE. So it would be the best to look for the energy loss effects by studying the 
production of hadrons in the nucleus fragmentation region. Quadratic dependence on L will 
be manifest in the A-dependence of the number of knocked out nucleons, as well as in the 
fluctuations of the number of emitted soft protons and neutrons in the case of heavy nuclei. 
Current HERA detectors have good acceptance for such nucleons. 

Broadening of the transverse momentum spectrum may be more easy to access. The trans- 
verse momentum of a parton increases as a result of multiple collisions. Since the momentum 
transfers in the subsequent collisions are independent, the increase of the (average) transverse 
momentum squared is proportional to L and hence to the number of secondary collisions. One 
expects for the quark 

A P ^0.2GeV^, (34) 

while for gluons broadening is about a factor of 2 larger. 

A word of caution is necessary here. For distances typical even for heavy nuclei the average 
momentum transfer is rather small so application of perturbative QCD may be difficult to 
justify. Also it is not clear whether it is safe to interpolate from fixed target energies to collider 
energies assuming that the momentum transfer in individual collisions is energy independent. 

To study these effects at HERA in a clean way one needs to consider processes dominated 
by relatively large x > 0.1. They include processes of dijet production similar to those studied 
at FNAL An advantage of HERA is that it would be possible to use information about 
the decay of the nucleus to check the correlation between transverse momentum broadening 
and the number of struck nucleons. 



3.7.4 Parton Propagation at Small x 

This region is more relevant for HERA but, unfortunately, theoretical estimates are difficult and 
uncertain because the system traversing the nucleus is fairly complicated. Hence its interaction 
with nuclear matter not easy to evaluate. Perturbative QCD leads to the simple prediction 
that at large Q 2 and large incident energies, due to QCD factorization, the spectrum of leading 
hadrons |72| is given by, 

1 drr^ +A ^ h+x (r O 2 z) 

_ML^ /feln(n (35) 
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which does not depend on A. (Energy losses we discussed above do not change the z spectrum 
in this limit). However, it is far from clear if perturbative methods are applicable at all - 
even at large Q 2 pi. The incoming system can experience several soft interactions with the 



nucleons of the target nucleus. However, the AGK technique which can be used to relate nuclear 
shadowing with the A dependence of diffraction and the fluctuations of hadron production in 
the central region, does not allow any predictions for the A dependence of the spectrum of 
leading hadrons in the current fragentation region since these effects depend on the details of 
the virtual photon wave function. In view of these difficulties it is not possible to give unbiased 
quantitative predictions and we restrict ourselves to a discussion of qualitative expectations 
and possible interpretations of the future measurements. 

As we have already mentioned, at small x the virtual photon fluctuates into a qq pair a long 
time before it enters the nucleus. This has several consequences. First, it opens the possibility 
of coherent phenomena in which the nucleus participates as a whole which were discussed 
in a previous section. They should be carefully separated from the incoherent interactions 
we are concerned with. Second, there is enough time for this fluctuation to emit a large 
number of gluons and new qq pairs before it enters the nucleus. The cross-section of such a 
"dressed" fluctuation is generally fairly large, leading to substantial "shadowing" effects, as 
already discussed. Finally, the system produced in the first collision is by no means a single 
bare quark but rather a multiparton conglomerate - result of a quark-gluon cascade of length 
equal to the available rapidity, i.e. very long at small values of x. Studies of interactions in 
nuclear matter provide an opportunity to obtain information about this object. 

The cascade origin of the system travelling through the nucleus implies strong correlations 
between partons. Consequently, large fluctuations are expected in its physical properties and 
thus also in the observable nuclear effects. For example, the cross section for secondary in- 
teractions is expected to vary widely from event to event. This variation may or may not 
be correlated with the fluctuations of the parton multiplicity and energy distributions. For a 
clean interpretation of the results it is therefore crucial to obtain information on the number of 
secondary interactions inside the nucleus which is the most important parameter determining 
the strength of the interaction. Fortunately, as we have already explained, this information is 
accessible through measurements of the distribution of nucleons (protons and neutrons) emitted 
by the colliding nucleus. They seem thus crucial for the success of this investigation. 

The distribution of the number of collisions gives straightforward information on the fluc- 
tuations of the cross section for secondary interactions inside the nucleus. Its change as Q 2 
increases from to the deep inelastic region should give information on the extent to which the 
structure of parton bound states (i.e. hadrons) differs from that of the "dressed" qq pair. 

Similar remarks apply to measurements of the energy loss (which measures the energy 
distribution inside the parton system in question). Correlation between the observed energy 
loss and the number of collisions gives information on clustering phenomena inside the parton 
system. It will be interesting to look for possible effects of constituent quarks at low Q 2 



and to see when they disappear as Q increases. Fig. ^0] illustrates the possibility of such 
an investigation. It shows the ratio of the total energy loss of the incident virtual photon in 
the nucleus to that in collision with one nucleon, assuming that the photon fluctuates into K 
constituents which interact independently (losing energy) with the target. One sees that the 
curves corresponding to different values of K are substantially different. This seems to give a 
chance to see possible effects of constituent quarks (K = 2) at low Q 2 which should gradually 
disappear (K — > oo) as Q 2 increases. It should be remembered, however, that variations in 
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Porton energy loss in a nucleus compared to a free nucleon 
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Figure 20: Energy loss of the incident virtual photon in nuclei assuming that the photon 
fluctuates into K constituents interacting independently in nuclear matter with an interaction 
cross section of 20 mb. 



the cross-section of the photon fluctuation in question will smear out this clean effect. It is 
therefore necessary to rely on an independent estimate of the number of collisions as we already 
emphasized at the beginning of this section. In this case the ratio in question is given by 

-^ = K(1-(1- l/Kf) (36) 

(N is the number of collisions) and one can easily see that this is indeed a much better way to 
analyze the data, particularly at high N. 

It should be noted at this point that, in contrast to the situation at finite x, in the present 
case the energy loss is expected to be a finite fraction of total available energy, as explained 
above. Consequently, there seem to be no particular difficulties in this measurement. 

To summarize, one expects a dramatic increase of the energy loss in nuclear matter when 
one goes from finite to small values of x. If confirmed by future data this should allow to study 
details of the parton structure of the QCD cascade. Investigation of the transition region of 
x ~ 0.01 is also of great interest. HERA is very well suited for this task. 



3.8 Connection to Heavy-ion Collisions at High Energies 

The interplay between the physics which can be studied in high-energy eA collisions at HERA 
and that to be studied in the heavy ion physics was discussed at the dedicated workshop 
"Nuclei at HERA and Heavy Ion Physics" which was held at Brookhaven National Laboratory 
in 1995. It was concluded that the measurements of eA collisions at HERA can provide crucial 
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information necessary for unambigous interpretation of the heavy ion colllisions at RHIC and 
LHC for establishing whether a quark-gluon plasma is formed in these collisions. 

Three major links are 

• Nuclear gluon shadowing 

One needs xg A (x,Q 2 ) for x ~ 1(T 2 , Q 2 ~ 1 - lOGeV 2 and x ~ lO" 3 , Q 2 ~ lOGeV 2 to 
fix the initial conditions at RHIC and LHC respectively. This is especially important for the 
LHC since mini-jet production determines the initial conditions for y/s > lOOGeV A. The bulk 
of the particles produced at central rapidities in AA collisions at the LHC is expected to be 
generated due to this mechanism p5 |. Currently uncertainties in nuclear shadowing transform 



into at least a factor 2-4 differences in the final transverse energy flow |96 |. 
• Jet quenching 

Recent QCD studies 0] have demonstrated that the medium induced energy losses and 



Pt broadening of a high energy parton traversing a hot QCD medium are much larger than in 
the case of a cold medium. This provides a unique new set of global probes of the properties 
of the state formed during AA collisions [95|. To interpret unambiguously this effect it is 



necessary both to measure the nuclear gluon shadowing and to study the parton propagation in 
cold matter in DIS to confirm that the energy losses (pt-broadening) remain small at energies 
comparable to those to be studied at RHIC and LHC. 

• Testing of soft dynamics of interactions with nuclei 

Study of eA interactions at HERA in the same energy range as that to be studied 
in pA and AA collisions at RHIC (\/s ~ 200 GeV) will provide a unique testing ground for 
the modern models of interactions with nuclei which aim at describing on the same footing 
ep, eA, pp, pA, AA collisions |89| . It would allow to be established whether or not the same 



dynamics determines hadroproduction in eA collisions and in central AA collisions. 



4 Acknowledgements 

We thank all our colleagues in working group 8 for their invaluable assistance and for numerous 
discussions. AB acknowledges support by the KBN grant No 2 P03B 083 08 and by a PECO 
grant from the EEC Programme "Human Capital and Mobility", Network "Physics at High 
Eneregy Colliders", Contract No ERBICIPDCT940613. MS thanks the U.S. Department of 
Energy for financial support under grant number DE-FG02-93ER-40771. 



References 

[1] V. N. Gribov, Sov. Phys. JETP 30 709 (1969). 

[2] J. D. Bjorken in Proceedings of the International Symposium on Electron and Photon 
Interactions at High Energies, p. 281-297, Cornell (1971). 



[3] R. Baier, Yu.L. Dokshitzer, A.H. Mueller, S. Peigne and D. Schiff, [hep-ph 9608321 and 
these proceedings. 



36 



[4] J. Chwastowski and M.W. Krasny, "What can we gain by detecting nuclear fragments in 
electron-nucleus collisions at HERA?" , these proceedings. 

[5] K. Kurek, "QED Radiative Processes in Electron-Heavy Ion Collisions at HERA", these 
proceedings. 

[6] I. Akushevich and H. Spiesberger, "Radiative Corrections to Deep Inelastic Scattering on 
Heavy Nuclei at HERA" , these proceedings. 

[7] NMC, M. Arneodo et al., Nucl. Phys. B 441 (1995) 3. 

[8] NMC, M. Arneodo et al, Nucl. Phys. B 441 (1995) 12. 

[9] E665, M.R. Adams et al, Z. Phys. C 67 (1995) 403. 

[10] NMC, M. Arneodo et al., "The A dependence of the nuclear structure function ratios", 
submitted to Nucl. Phys. B. 

[11] NMC, M. Arneodo et al., "The Q 2 dependence of the structure function ratio F^/F^ and 
the difference R Sn — R c in deep inelastic muon scattering" , submitted to Nucl. Phys. B. 

[12] EMC, J.J. Aubert et al, Phys. Lett. B 123 (1983) 275. 

[13] EMC, J. Ashman et al, Phys. Lett. B 202 (1988) 603 and Z. Phys. C 57 (1993) 211. 

[14] EMC, M. Arneodo et al, Phys Lett B 211 (1988) 493 and Nucl. Phys. B 333 (1990) 1. 

[15] SLAC, A. Bodek et al, Phys. Rev. Lett. 50 (1983) 1431, ibid. 51 (1983) 534 and R. G . 
Arnold et al, ibid. 52 (1984) 727. 

[16] BCDMS, G. Bari et al, Phys. Lett. B 163 (1985) 282 and A.C. Benvenuti et al, ibid. B 
189 (1987) 483. 

[17] For a review of the experimental and theoretical situation of nuclear effects in DIS see e.g. 
L.L. Frankfurt and M.I. Strikman, Phys. Rep. 160 (1988) 235; 
T. Sloan, G. Smadja and R. Voss, Phys. Rep. 162 (1988) 45; 

R.J.M. Covolan and E. Predazzi, in "Problems of Fundamental Modern Physics", Editors 
R. Cherubini, P. Dalpiaz and B. Minetti, World Scientific, Singapore (1991) p. 85; 
M. Arneodo, Phys. Rep. 240 (1994) 301. 

[18] J. Kwiecinski,A.D. Martin and P.J. Sutton, Phys. Rev. D52 (1995). 

[19] A.H.Mueller and J-W.Qiu, Nucl. Phys. B268 (1986) 427. 

[20] A.D. Martin, R.G. Roberts and W.J. Stirling RAL Report TR 96-037 and Int. J. of Mod 
Phys A10 (1995) 2885. 

[21] HI, S. Aid et al, Nucl. Phys. B449 (1995) 3. 

[22] D. Graudenz, CERN-TH-95-121, April 1995. Talk given at Workshop on Deep Inelastic 
Scattering and QCD (DIS 95), Paris, France, 24-28 April 1995. 

[23] EMC, J.J. Aubert et al., Nucl. Phys. B 213 (1983) 1 and J. Ashman et al, Z. Phys. C 
56 (1992) 21. 

37 



[24] NMC, D. Allasia et al, Phys. Lett. B 258 (1991) 493. 
[25] NMC, P. Amaudruz et al., Nucl. Phys. B 371 (1992) 553. 
[26] HI, S. Aid et al, Nucl. Phys. B 472 (1996) 3. 

[27] ZEUS, M. Derrick et al., paper pa-02 47 submitted to the XXVIII International Conference 
on High Energy Physics Warsaw, July 25-31, 1996. 

[28] HI, S. Aid et al., paper pa-02 73 submitted to the XXVIII International Conference on 
High Energy Physics Warsaw, July 25-31, 1996. 

[29] T. Weiler, Phys. Rev. Lett. 44 (1980) 304. 

[30] E.L. Berger and D. Jones, Phys. Rev. D 23 (1981) 1521. 

[31] R. Baier and R. Riickl, Nucl. Phys. B 201 (1981). 

[32] A.D. Martin, C.-K. Ng and W.J. Stirling, Phys. Lett. B 191 (1987) 200. 

[33] M. Kramer et al, Phys. Lett. B 348 (1995) 657 and Nucl. Phys. B 459 (1996) 3. 

[34] ZEUS,M. Derrick et al. Phys. Lett. B315 (1993) 481. 

[35] HI, T. Ahmed et al. Nucl. Phys. B429 (1994) 477. 

[36] HI, T.Ahmed et al., Phys.Lett. B348 (1995) 681. 

[37] ZEUS, M. Derrick et al. Z. Phys. C68 (1995) 569. 

[38] H. Jung, "Hard Diffractive Scattering in High Energy ep Collisions and the Monte Carlo 
Generator RAPGAP", DESY preprint DESY 93-182 (1993), to appear in Comput. Phys. 
Commun. 

[39] W. Buchmuller and A. Hebecker, Phys.Lett. B355 (1995) 573. 
[40] L.Frankfurt and M.Strikman, Phys.Lett. B382 (1996) 6. 
[41] M.G. Ryskin, Z. Phys. C 57 (1993) 89. 

[42] S.J.Brodsky, L.Frankfurt, J.F.Gunion, A.H.Mueller and M.Strikman, Phys. Rev. D50 

(1994) 3134. 

[43] J. Nemchik et al., Phys.Lett. B341 (1994) 228. 

[44] L. Frankfurt, W. Koepf and M. Strikman, Tel-Aviv University preprint TAUP-2290-95 

(1995) ; Phys.Rev. D54 (1996) 3194. 

[45] M.G. Ryskin et al., "Diffractive J/ip photoproduction as a probe of the gluon density", 
Rutherford Lab. preprint RAL-TR-95-065, |hep-ph/9511228| (1995). 

[46] M. Arneodo, L. Lamberti and M. G. Ryskin, DESY report DESY 96-149, to appear in 
Comp. Phys. Comm. 

[47] ZEUS Collab., M. Derrick et al, Phys. Lett. B 356 (1995) 601. 

38 



[48] EMC, J. Ashman et al, Z.Phys. C52 (1991) 1. 
[49] E665, M.R.Adams et al., Z.Phys. C61 (1994) 179. 
[50] N.Pavel, these proceedings. 

[51] M. Krawczyk and B. B. Levtchenko, these proceedings. 

[52] B. Blattel, G. Baym, L.L. Frankfurt and M. Strikman, Phys. Rev. Lett. 71 (1993) 896. 
[53] L. L. Frankfurt and M. Strikman, Phys. Rep. 160 (1988) 235. 
[54] L.M Lerran and R.Venugopalan, Phys.Rev. D50 (1994) 225. 

[55] V.Abramovskii, V. N. Gribov and O. V. Kancheli, Sov. J. Nucl. Phys. 18, 308 (1974). 
[56] L. L. Frankfurt and M. I. Strikman, Nucl. Phys. B316 (1989) 340. 

[57] B.Z.Kopeliovich and B.Povh, Phys.Lett. B367 (1996) 329 and contribution to these pro- 
ceedings. 

[58] B.K. Jennings and G.A. Miller Phys. Rev. C49 (1994) 2637. 
[59] HI, S.Aid et al., Z. Phys. C69 (1995) 27. 

[60] T.Gousset and H.J. Pirner, Phys. Lett. B375 (1996) 354 and contribution to these proceed- 
ings. 

[61] L.Frankfurt, M.Strikman and S. Liuti, Phys. Rev. Lett. 65 (1990) 1725. 
[62] J.W.Qiu Nucl.Phys. B291 (1987) 746. 

[63] L. Frankfurt, M. Strikman and S. Liuti, Proceedings of the Conference on Particles and 
Nuclei XIII, Perugia, Italy (1993) page 342, edited by A. Pascolini (published by World 
Scientific). 

[64] K.J. Eskola, Nucl.Phys.B400 (1993) 240. 

[65] K.J. Eskola, Jian-wei Qiu, Xin-Nian Wang, Phys.Rev.Lett.72 (1994) 36. 
[66] N.N.Nikolaev and V.G.Zakharov, Z.Phys. C49 (1991) 607. 

[67] A.H.Mueller, In the proceedings of the BNL workshop on Nuclei at HERA and Heavy Ion 
Physics, BNL-62634 (1995). 

[68] A.H.Mueller and H.Navelet, Nucl.Phys. B282 (1987) 727. 

[69] A.Donnachie and PLandshoff, Phys.Lett B296 (1992) 227. 

[70] F. E. Low, Phys. Rev. D12 (1975) 163. 

[71] S.J. Brodsky in Proceedings of the Thirteenth International Symposium on Multiparticle 
Dynamics, ed. W. Kittel, W. Metzger and A. Stergiou (World Scientific, Singapore, 1982) 
page 963. 



39 



[72] A.H. Mueller in Proceedings of the Seventeenth Rencontre de Moriond, Moriond, 1982 ed. 
J. Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, France, 1982) Vol. I, page 13. 

[73] E665,M.R. Adams et al, Phys.Rev.Lett. 74 (1995) 1525. 

[74] L. Frankfurt, G.A. Miller and M. Strikman, Phys. Lett. B304 (1993) 1. 

[75] J.Bartels,H.Lotter and M. Wusthoff, Phys.Lett. B379 (1996) 239. 

[76] J.Bartels, C.Ewerz, H.Lotter, M.Musthoff, DESY 96-085. 

[77] H.Abramowicz, L.Frankfurt and M.Strikman, DESY-95-047; SLAC Summer Inst. 1994:539- 
574. 

[78] L. Frankfurt, W.Koepf, M.Sagrsyan, M.Strikman, contribution to these proceedings. 

[79] ZEUS, M.Derrick et al., Phys.Lett. B369 (1996) 55. 

[80] J.D.Bjorken, Phys.Rev. D47 (1992) 101. 

[81] A.H.Mueller and W.-K.Tang, Phys.Lett. B284 (1992) 123. 

[82] V.Del Duca and W.-K.Tang, Phys.Lett. B312 (1993) 225. 

[83] DO, S.Abachi et al., Phys.Rev.Lett. 72 (1994) 2332; FERMILAB-PUB-95-302-E(1995). 

[84] CDF, S.Abe et al,Phys.Rev.Lett. 74 (1995) 85. 

[85] L.Frankfurt and M.Strikman, contribution to these proceedings. 

[86] B.Andersson, I.Otterlund and E.Stenlund, Phys.Lett. B73 (1978) 343. 
J.Babecki and G.Nowak, Acta Phys. Polonica B9 (1978) 401. 

[87] see e.g. J.Babecki et al., Acta Phys.Polonica B9 (1978) 495. 

[88] E665, M.R.Adams et al Phys.Rev.Lett. 74 (1995)5198. 

[89] K. Geiger, Contribution to these proceedings; J.Ellis, K. Geiger and H. Kowalski, to be 
published in Phys.Rev.D, e-Print Archive: [hep-ph / 9605425 . 

[90] M.Strikman, M.Tverskoy, M.Zhalov, these proceedings. 

[91] D.M. Aide et al., Phys.Rev.Lett. 66 (1991) 2285 and references therein. 

[92] T. Fields and M.D. Corcoran, in the EPS Conference proceedings, Marseille, France, July 
22-28, 1993; Phys.Rev.Lett. 70,143,1993; 
R.C. Moore, et al, Phys.Lett. B244 (1990) 347 ; 
M.D. Corcoran, et al., Phys.Lett. B259 (1991) 209. 

[93] M. Luo, J. Qiu and G. Sterman, Phys.Rev.D49(1994) 4493. 

[94] J.Bjorken, SLAC-PUB-7096, January 1996. 

[95] X-N.Wang and M.Gyulassy, Phys.Rev.Lett. 68 (1992) 1480. 

[96] M.Gyulassy, Proc. 'Nuclei at HERA and Heavy Ion Physics', BNL-62634, 1995. 

[97] R. Baier, Yu.L. Dokshitzer, A.H. Mueller, S. Peigne and D. Schiff |EeFph 9607355 . 



40 



